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Abstract 

We show that a relativistic gas may be at "global" equilibrium in the 
expanding universe for any equation of state < p < p/3, provided that 
the gas particles move under the influence of a self-interacting, effective 
one-particle force in between elastic binary collisions. In the force-free limit 
we recover the equilibrium conditions for ultrarelativistic matter which 
imply the existence of a conformal timelike Killing vector. 
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1 Introduction 



The Lie derivative £^^gab of the metric gab with respect to may be used to 
characterize spacetime symmetries. The vector is a Killing vector (KV) if 
£i,9ab = is satisfied. In fluid spacetimes the vector ^* = u^/T, where 
is the four-velocity of the fluid, normalized according to u^ui = —1, and T is 
its equilibrium temperature, is of particular interest. In case n*/T is a KV, 
the corresponding spacetime is stationary. In a cosmological context one may 
consider the possibility of u*/T to be a conformal Killing vector (CKV): 

£ua9ik=(^] + ( ^ ) = '^(t>9^k ■ (1) 



It is easy to realize that the latter relation implies 
as well as 



i k n T @ d 

-uu £^^gik = — = -— = , (3) 
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where T = T^au"" and G = u?j is the fluid expansion. The quantity a is a 
length scale which in the special case of a homogeneous and isotropic universe 
coincides with the scale factor of the Robertson Walker metric. Equation (^) 
implies T oc which describes the temperature behaviour of an ultrarela- 
tivistic fluid (radiation) in the expanding universe. The ultrarelativistic limit 
is clearly singled out by the CKV property, although no equation of state was 
used to obtain the temperature behaviour (^). It is well known too, that for a 
relativistic Boltzmann gas the CKV property of u^'/T follows as an equilibrium 
condition ("global" equilibrium) from Boltzmann's equation in the correspond- 
ing limit of massless particles (radiation) |2|, ^ For any other equation of 
state the "global" equilibrium condition is £ua/T9ik = 0, a requirement which 
cannot be realized in an expanding universe. The ultrarelativistic equation of 
state is a singular case both macroscopically and microscopically. This cir- 
cumstance may give rise to the question as to whether there exist symmetries 
which are compatible with a nonvanishing expansion and which prefer different 
equations of state in a similar sense in which the CKV symmetry prefers an 
ultrarelativistic equation of state. It is this question which we will address and 
answer affirmatively in the present paper. 

Based on the observation that the temperature law (^), which is implied 
by the CKV property (g) of m/T, may independently be derived from the 
thermodynamics of an ultrarelativistic fluid, we find, in section 2, an expression 
for £ua/Tgik which fits the temperature law for a fluid with arbitrary equations 
of state. Section 3 shows how this expression follows from the Boltzmann 
equation as an equilibrium condition for the one-particle distribution function 
for a relativistic gas in a force field. In section 4 we consider the microscopic 
particle motion in this force field and confirm explicitly that the equilibrium 
distribution of the gas particles is maintained as the universe expands. Units 
have been chosen so that c = A;^ = 1. 

2 Spacetime symmetries 

In a first step we recall how the temperature law (^), which we obtained here 
as a consequence of the CKV property (g), is alternatively and independently 
derived on a thermodynamical basis. Let us consider a perfect fluid with the 
energy-momentum tensor 

T^'' = pu'u^ + pK'^ , (4) 

where p is the energy density measured by a comoving observer, p is the equi- 
librium pressure, and h^^ is the projection tensor h^^ = g^^ + u^u^ . Local 
energy-momentum conservation T*^ = implies 

-UiT% = p + Q{p + p) = ^ (5) 

and 

himT% = {p + p)u^ + VmP = . (6) 
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The number N of particles in a comoving volume is = na^, where n is 
the particle number density. The corresponding particle number flow vector is 
N"" = nu"" and the particle number conservation law may be written as 

AT.^ = n + en = . (7) 

Given equations of state in the general form 

p = p{n,T) , p = p{n,T) , (8) 

differentiation of the latter relation and using the balances (|5|) and (0) yields 

^ = -e^ (9) 

dp _ (dp/dT)^ 



where the abbreviation 



dp {dp/dT)^ 



as well as the general relation 

dp\ p + P T / dp\ 
dn J T n n \ dT J „ 

have been used. For an ultrarelativisic equation of state dp/ dp = 1/3 we recover 
the temperature law (^). It is only this specific equation of state for which the 
expressions (|^ and (|3|) are compatible. 

It is now tempting, however, to combine the temperature behaviour (^) with 
the Lie derivative of the metric also for a general equation of state. It is easy 
to realize that for a choice 



£ ua 9ik = 2</> 



gik - ( 3^ - 1 ) UiUk 



(10) 



where is again given by (^), the temperature law (^) is reproduced for any 
equation of state. In the ultrarelativistic limit dp/ dp = 1/3 we recover the CKV 
case (^). The expression ( p!o| ) for the Lie derivative makes the fluid spacetime 
symmetry depend on the fluid equation of state. The CKV property ^ is no 
longer a singular case. 

For a gas, the KV and CKV conditions are known to follow from Boltz- 
mann's equation as equilibrium conditions for the one-particle distribution func- 
tion [|l|, ^, 1^, 1^. In the following section we will clarify the conditions under 



which the general expression (10) is related to equilibrium properties of the gas. 
The essential point will be that in order to obtain ( p!o| ) as equilibrium condition 
for the one-paricle distribution function, the gas particles will have to move 
in a suitable force field. The ultrarelativistic case, corresponding to the CKV 
property, will then follow in the limit of a vanishing force. 



3 



3 Kinetic theory in a force field 



Let the particles of a gas in between elastic binary collisions move under the 
influence of a four force = -F* {x,p) according to the equations of motion 

<^=p-, 2^=mF'. (11) 

d7 d7 

where 7 is a parameter along their worldline which for massive particles may 
be related to the proper time r by 7 = r/m. Since the particle four-momenta 
are normalized according to p^pi = —m?, the force -F* has to satisfy the relation 
PiF'' = 0. The one-particle distribution function / = f{x,p) then obeys the 
Boltzmann equation (see, e.g., ^) 

p7,.-rLpV0 + mF^0 = C[/] . (12) 

C[f] is Boltzmann 's collision integral. The effective one-particle force is 
supposed to model different kinds of interactions in a simple manner |^, |9|. 
Its specific structure in the present context will be determined below. 

The particle number flow 4-vector A^* and the energy momentum tensor T**^ 
are defined in a standard way (see, e.g., ^) as 

N' = J dPp'fix,p) , r' = J dPp'p'^fix.p) . (13) 



The integrals in the definitions (|T3|) and in the following are integrals over the 
entire mass shell p^pi = —m?. The entropy flow vector 5" is given by |^, ^] 

5'^ = - jp^lflnf-ndP , (14) 

where we have restricted ourselves to the case of classical Maxwell-Boltzmann 
particles. 

Using well-known general relations (see, e.g., [||) we find 



and 



= Jp''[c[f]-mF^^^dP , (15) 

S^, = -j\nf{c[f]- dP . (16) 

Since we are interested in the equilibrium properties of the gas we will restrict 
ourselves to collisional equilibrium from now on. Under this condition In / in 
Eq. ( [l6| ) is a linear combination of the collision invariants 1 and p"" and the 
collision integral C [/] vanishes. The corresponding equilibrium distribution 
function becomes (see, e.g., Q) 

/0(x,p) = exp[a + /?,p'^] , (17) 
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where a = a [x) and (ia (x) is timelike. Inserting the equihbrium distribution 



function (|17| ) into the Boltzmann equation (12) one obtains 

P"«,a + /3(a;6)Py = -m^F^ . (18) 

For a vanishing force the latter condition reduces to the "global" equilibrium 
condition of standard relativistic kinetic theory, i.e., to a = const and either to 
the KV condition for m > 0, or to the CKV condition (|l|) for m = 0. Use of 
the equilibrium distribution function (|^) in the balances ( |T5|) yields 



N^, = -mPi J F'fdP , 



T'^^ = -m(5i j p^F'fdP . (19) 

For the entropy production density (|l6|) we find 

= + PaP"] F'fdP = -aN^^ - paT% . (20) 

With / replaced by f in the definitions (|l|) and (0), A^'', T''^ and S" may 
be split with respect to the unique 4-velocity u"" according to 

N" = nu" , T"^ = pu'^u'' + ph"'' , S" = nsu" , (21) 

where n", h""^, n, p, and p may be identified with the corresponding quanti- 
ties of the previous section, s is the entropy per particle. The exact integral 
expressions for n, p, p are given by (see, e.g., Q), 



p AT:m?T ^ / m 



n = ^ = ^:r^K2 ( - ) exp [a] , (22) 



and 



= = ;; = »lfjtfj + 3T = .n^-T. (23) 

The quantities Kn are modified Bessel functions of the second kind The 
entropy per particle s is 

s=P±?-a. (24) 



In order to evaluate the equilibrium condition (|18|) , we expand the quantity 
/?jF* in a power series in p"" . According to the structure of the left-hand side of 
(|l^), it will be sufficient to consider terms up to second order in the momentum: 

AF* = AF>'^ + AF|,/y • (25) 

The quantities F^ and F^; are spacetime functions to be determined below, i.e., 
they do not depend on p"" . We emphasize that it is only the projection /JjF' of 
F* which is expanded in a power series, not the force itself. The orthogonality 
of the latter to p^ may always be guaranteed [§]. 
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With the expansion (p5| ) the equihbrium condition (|T8|) sphts into 

-mPiFi (26) 



a., 



and 

With the identification (3i = m/T it is obvious that the choice 



ab 



e_ 

ST 



9ab 



(27) 



(28) 



just reproduces (10). The symmetry requirement ( [7^ is recovered as an equi- 
librium condition for the one-particle distribution function for gas particles in 
a specific force field. For dp/ dp = 1/3 the condition (^) with (28) reduces to 
the CKV condition @. 

The condition (p6|) may be written as 



- mpiFl = -aua + Vaa , (29) 
where V^a = /i^a,;,. Applying the ansatz (p^) in the balances (p^), we find 



and 



rpak 



-mA {fI 



(30) 



(31) 



where M"^^ = J dP/° 

pO-pkpl jg third moment of the equilibrium distribution 
function. Use of the expressions (^) and ( p9| ) allows us to write the particle 
number balance (^) as 



na + nQ 



1 



p dp 
nT dp 



(32) 



With the identification a = fi/T, the quantity a may be calculated from the 
Gibbs-Duhem relation 



dT 



dp={p + p) — + nTd [ - 



With p = nT and the temperature law (^) we obtain 

= -G 



a = \ — 



1 _ P 
nT dp 



(33) 



(34) 



Inserting this result into (32) we recover formula (|^) i.e., the particle number 
is indeed conserved. Analogously, we find from Eq. (pl[). 



rpak 
Ual -k 



-ap + (f)Ua 



9kl 



'dp ~ ' 



M' 



akl 



Obviously, one has 



(35) 



(36) 
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The expression UaUkUiM"'^^ in the balance ( ^5| ) may be evaluated by use of the 
relations 



a rO 



(the derivative has to be taken for a = const), 

^ = -T{p + p) and ^ = -3T {p + p) - pT 



(37) 



which follow with the help of the properties 



dz \ z" / z" 
of the functions Kn H. The result is 



and Kn+i = 2n h Kn-i 

z 



UaUkUiM 



akl 



-3T{p + p)-pT 



The energy balance (|3q ) then becomes 



dp 
dp 



1 + 5 



p + p f P + P 



nT 



nT 



nTQ 



Since the expression in the bracket of the last relation is H 



t) 



nT 



P + P 
nT 



dp 

dp 



(38) 



(39) 



(40) 



(41) 



(42) 



one recovers the local energy conservation (|5|). Similarly one shows that the 
momentum balance (m) is fulfilled identically, which implies 



VaO; 



Ua + 



VaT\ p + p 

nT 



. 



Thus we have clarified that the self- interacting force term ( psD with (28) and 
(p9|) both implies the symmetry condition ( [To| ) and at the same time respects 
the conservation laws for particle number and energy momentum. As a conse- 
quence, the entropy production density (pO|) vanishes also. We emphasize that 
the mentioned conservation properties single out a very specific self-interaction. 
For a general force term of of the type (|2^) with a coefficient different from (^) 
none of the moments of the distribution function (|l7|) will be conserved. This 
has been explicitly demonstrated for two cases, the first of which is a confor- 
mal symmetry (0) for massive particles |^ 0, While such a configuration 
does not exist for a conserved particle number, it turned out to be possible if 
these particles were created at a specific rate. The second case was based on 
a modification of the conformal symmetry (||) according to £ua/T9ik = hik 
The corresponding self-interaction implies particle production at a rate which 
coincides with the expansion rate of the universe.. The source terms in the 
balances for the second moments in each of these cases were mapped onto bulk 
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pressures of effective energy momentum tensors which are conserved. Although 
self-interactions of this kind are entropy producing, the corresponding states 
are "generalized" equilibrium states ||^ since the gas particles are still charac- 
terized by the equilibrium distribution function (17). The case dealt with here 
is singled out by the fact that the "generalized" equilibrium coincides with the 
conventional "global" equilibrium with conserved moments of the distribution 
function. 



4 Microscopic particle motion 



Knowledge of the coefficients ( pq ) and ( [29| ) in ( pq ) allows us to study the motion 
of the gas particles explicitly. Contracting the equation of motion in (|ll]) with 
the macroscopic four-velocity results in 



^ ' = UiF' + —Ui-kPp 
m 



dr 



(43) 



where we have used that 



dr ' m 



With the well-known decomposition of the covariant derivative of the four- 

(44) 



velocity [|o|, |T|], 



G 



where aab = KK'^ {c;d) - KbQ/^ and ujab = KH'^lc;d\, we find generally 



D (uip*) 



UiF' + —Qhikp'p'' + —(Tikp'p 
6m m 



1 



dr 



1 



m 



iiiUkp'p'^ 



(45) 



In the following we investigate this equation for the homogeneous case Ui = 0. 
(Notice that aab = Ua + ^aT/T = are consequences of (p^)). For UjF* we 
have under such conditions 



UiF' = - 
m 



e e 



aE + m^— + —E^ 3-f- - 1 



3T 3T 



dp 



dp 



(46) 



where E = —Uip^. Equation ( |45| ) simplifies to 

— = -UiF' [E^ - 

dr 3m V 



dr ' 3m 

With dr = dt [m/E) and dE/dt = E, the last equation is equivalent to 



(47) 



[E^ - m?- 



,2V 



+ 



2m 



TUiF' 



(48) 



^2 _ j^2 g2 J^2 _ 

In the force- free case F^ = 0, i.e., for geodesic particle motion, we find 

(r = 0) , (49) 



TT^2 2 -2 

E — m (X a 
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implying the expected behavior E (x a"^ for massless particles (photons), while 
the nonrelativistic energy e = E — m with e ^ m of massive particles decays 
as e oc a~^. 

In the presence of the force term (p6|) , use of (^) and dH) allows us to obtain 

= A + (ii; - /i) ^ , (50) 

equivalent to 

— = f = a + PaP = const . (51) 

This demonstrates explicitly that the equilibrium distribution (|l^) is indeed 
maintained in the expanding universe for arbitrary equations of state. 



5 Conclusion 

We have clarified the conditions under which the particles of an expanding 
gaseous fluid are governed by an one-particle equilibrium distribution func- 
tion which satisfies the relativistic Boltzmann equation. In general, such kind 
of equilibrium is only possible if there exist suitable forces inside the system 
which make the particles move on nongeodesic trajectories in between equi- 
librium establishing elastic binary collisions described by Boltzmann's collision 
integral. These forces depend on the microscopic particle momenta as well as 
on macroscopic fluid variables including the equation of state. They represent 
a (nongravitational) self-interaction of the gas which is consistent with the con- 
servation of particle number and energy momentum. In other words, the gas as 
a whole must mobilize specific internal forces to maintain an equilibrium dis- 
tribution of the particles as the universe expands. The familiar ultrarelativistic 
case in which Ui/T is a CKV corresponds to the force free limit. 
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